INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE 



Henry' law and gas phase disappearance 

Jerome Jaffre — Amel Sboui 



N° 6891 

Mars 2009 

Theme NUM 



INSTITUT NATIONAL 
DE RECHERCHE 
EN INFORMATIQUE 
ET EN AUTOMATIQUE 




INKIA 



centre de recherche 

PARIS - ROCQUENCOURT 



Henry' law and gas phase disappearanc^ 

Jerome Jaffre0 , Amel SbouiH 

Theme NUM — Systemes numeriques 
Equipe-Projet Estime 

Rapport dc recherche n° 6891 — Mars 2009 — [5] pages 



Abstract: For a two-phase (hquid-gas) two-component (water- hydrogen) sys- 
tem we discuss the formulation of the possible dissolution of hydrogen in the 
liquid phase. The problem is formulated as a set of nonlinear partial differential 
equations with complementary constraints and we show how Henry's law fits in 
a phase diagram. 
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La loi de Henry et la disparition de la phase 

gazeuse 

Resume : Pour un systeme diphasique (liquide-gaz) a deux composants (eau- 
hydrogcnc) on considere la formulation de la dissolution eventuelle de I'hydrogcne 
dans la phase aqueuse. On montrc comment la loi de Henry entrc dans un dia- 
gramme de phase. Le probleme est alors formule comme un ensemble d' equations 
aux derivees partielles non-lineaires avec des contraintes de comp[lementarite. 

Mots-cles : Milieu poreux, ecoulement diphasique, dissolution, stockage pro- 
fond de dechets nucleaires 
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1 Introduction 

Recently the Couplex-Gas benchmark (http: / /www. gdrmomas.org/ex_qualifications. html) 
was proposed by Andra and MoMAS in order to improve the simulation of the 
migration of hydrogen produced by the corrosion of nuclear waste packages in 
an underground storage. This is a system of two-phase (liquid-gas) flow with 
two components (hydrogen- water) . This is similar to a black-oil model but 
the benchmark was proposing to use in Henry's law to model the dissolution 
of hydrogen. This benchmark generated some interest and engineers encoun- 
tered difficulties in handling the appearance and disappearance of the phases. 
Recently two other papers [5], [I] adressed the formulation of this problem. 

In this paper, for the sake of the simplicity of exposition, we will assume that 
only hydrogen is in the gas phase (water does not vaporize). We discuss how to 
extend Henry's law to the case where the gas phase disappears. This is done by 
showing how Henry's law fits in a phase diagram. The resulting formulation is a 
set of partial differential equations with complementarity constraints. A similar 
formulation is also presented in [5] and it also can be viewed as a variation of 
the formulation of the black oil model presented in 

2 Mathematical formulation 

We first describe the equations and some simplifying hypotheses. Since we 
consider a problem where the gas phase can disappear while the liquid phase 
cannot disappear we will consider a formulation using the saturation and the 
pressure of the liquid phase as the two main unknowns. 

2.1 Fluid phases 

£ and g are the respective indices for the liquid phase and the gas phase. Darcy's 
law reads 

q, = -K{x)h{s,){yp,- p.gVz), i^£,g, (2.1) 

where K is the absolute permeability. Assuming that the phases occupy the 
whole pore space, the phase saturations Si,i = £, g satisfy 

< Si < 1, S( + Sg — 1. 

The mobilities are denoted by ki{si) = — — with fc^i the relative permeability, 

Pi 

and Pi the viscosity of the phase i = £,g, assumed to be constant. The mobility 
ki is an increasing function of Si such that fci(O) = O^i = i, g. 

The phase pressures are related through the capillary pressure law 

Pc{si)=Pg~Pe>0. (2.2) 

assuming that the gas phase is the nonwetting phase. The capillary pressure is 
a decreasing function of the saturation S£ such that Pc(l) — 0. 

In the following we will choose Si,pi as the main variables since we assume 
that the liquid phase cannot disappear for the problem under consideration. 
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2.2 Fluid components 

We consider two components, water and hydrogen denoted by the indices j = 
w,h. The mass density of the phase i is 

where and p^j ^^re the mass concentrations of water and of hydrogen in the 
phase i,i = i, g. The mass fractions are 

Pi Pi 

Obviously we have 

We assume that the hquid phase may contain both components, while the gas 
phase contains only hydrogen, that is the water does not vaporize. In this 
situation we have 

Pa 

A third main unknown will be Xh- 

2.3 Conservation of mass 

We introduce the molecular diffusion fluxes jf,i = uu^h for the diffusion of 
component i in the liquid phase. They must satisfy + = and we have 

ji = ~cf>sePiDiVxl (2.3) 

The coefficients Df^ is a molecular diffusion coefficient. 

Conservation of mass applied to each component, water and hydrogen, gives 

d 

(t'^XsiPiXw) + div(p^xLqf + ji) = Qw, 
Q (2.4) 

(t^Q^isePexi + SgPg) + Ah7{pixi^i + Pg^g + ii) ^ Qh- 

We assume also that the gas is slightly compressible, that is pg = CgPg with Cg 
the compressibility constant, and that the liquid phase is incompressible, that 
is p^ is constant. 



3 Phase equilibrium 

3.1 Phase diagram 

A phase diagram is a figure like that shown in Figure [3Tl where a curve f{P, C) — 
separates two zones, a first zone defined by /(P, C) > where the systeme is 
monophasic - liquid or gas - and a second zone defined by /(P, C) < where 
the system is diphasic - liquid + gas. It gives the repartition of two components, 
like hydrogen and water, into the two phases, liquid and gas. 
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For the example we are concerned with, that is the eventual disappearance 
of the gas phase, we are interested in the left part of the phase diagram. There 
we have cither s; = 1 and /(P, C) > in the liquid zone, or s; < 1 and 
/(P, C) < in the diphasic zone. This can be written in terms of the following 
set of complementarity constraints 

(1 - s,)f{P, C) = 0, 1 - > 0, /(P, C) > 0. 

The variable C is the total concentration of hydrogen, that is in the diphasic 
zone C = + pg (since we assumed that only hydrogen can be present in the 
gas phase), in the liquid zone liquide C — pf^ and in the gas zone C = pg since 
we assumed that the water cannot vaporize. 

The variable P is called the " gas pressure" , P = in the gas and diphasic 
zone. To extend this definition to the liquid zone we express Pg in terms of the 
capillary pressure, pg = pi +pc{s{), and define P = pt +Pc(sf)- Thus in the 
liquid zone we have si = l,pc{si) = et P = p£. The definitions of the variables 
P and C are shown in Table [TJ 





liquid 


liquid + gas 


gas 




S£ = 1 


< se <1 


Si = 


p = 


Pi 


Pg =Pe+ Pc{si) 


Pa 


c = 


Pi 


pi + Pa 


Pa 


f{P,C) 


> 


< 


> 



Table 1: Definition of the variables P et C in the phase diagram. 



RR n° 6891 



6 



Jerome Jaffre and Amel Sboui 



3.2 Extension of Henry's law 

In the following we focus on the interface between the liquid zone and the 
diphasic zone. In the presence of the gas phase Henry's law reads 

Hp, = pi (3.1) 

with H the Henry constant. We can integrate this relation into a phase diagram 
formulation which includes the situation with no gas phase. This formulation 
includes a set of complimentarity constraints. 

Henry's law is defined only in the diphasic zone. In this zone we have 
HP — C = Hpg — {pi + Pg) = —pg < 0, and when we move to the curve 
separating the liquid and diphasic zones we have HP — C = pg — > 0. Above this 
curve is the liquid zone HP — C = Hpg — pi > 0. Thus using Henry's law is 
actually replacing the part of the curve f{P, C) ~ Q which separates the liquid 
and diphasic zones by the straight line HP — C = 0. This is shown in Figure [372] 
where only the left part of the phase diagram is shown with the line separating 
the liquid zone and the diphasic zone. 

Indeed, either the gas phase exists, 1 — > 0, Henry's law applies and 
H{pi +Pc{si)) — pi = 0, or the gas phase does not exist, sg = l,Pc{se) = and 
Hpi — pi > which says that for a given pressure pi the concentration is too 
small for the hydrogen component to be partly gaseous, or conversely for a given 
concentration pi the pressure p^ is too large for the hydrogen component to be 
partly gaseous. These cases can be written as the complementarity constraints 

(l-s,)(i/(p,+Pc(s^))-pl) =0, l-s,>0, Hipe+Pc{se))~pi>0. i^-^) 



P h 




C 

Figure 3.2: Part of the phase diagram separating the liquid and diphasic zones 
when Henry's law is used. 
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4 Formulation with complementarity constraints 

Finally we end up with a system of nonlinear partial equations - conservation 
equations and Darcy laws - with nonlinear complementarity constraints de- 
scribing the transfer of hydrogen between the two phases, the unknowns being 

d 

(l)-^{sipe{l - xD) + div(pf(l - xi)(li~ji) = Qw 
d 

fp-QlisePtxi + (1 - se)Cg{pe+pc{se))) + 

div(p£xlq£ + Cg{pi + Pc{si))cig + jl) = Qh 

= -K{x)kfXsi){Vpe - ptgVz) 

Qg = -K{x)kg{l - st){V[pi +Pc[st)) - Cg{pi +Pc{si))gVz) 

{l-Si){H{pe+p,{si))- pixi) =0, 1-S£>0, H{pe+p,{se))~ pexi>0- 

(4.1) 

This formulation is valid whether the gas phase exists or not. A similar system 
was presented in [2]. Wc note that when the gas phase disappears, i.e. S( ~ I, 
the second equation in (|4.1|) becomes a standard transport equation of hydrogen 
as a solute in the liquid phase. 

Assume that we discretize problem (|4.ip with cell centered finite volumes 
and we denote by N, the number of degrees of freedom for S£,pe,Xh which is 
equal to the number of cells. We introduce 

X € M.^^ , the vector of unknowns for se,pi,Xh, 

H : M.^-^ — > R^^, the vector of the discrctizcd component conservation equations. 

. ]g3N _^ ^Yic discrete vector for I - se, 
g : _^ ^ ^Yic discrete vector for H{pi + p^isi,)) - ptxi- 

Then the problem can be written in compact form 

n{x) = 0, 

T{x) '^g{x) = 0, T{x) > 0, g{x) > 0. 

Building a solver to solve properly such a problem is our current task. This 
task is undertaken in the field of numerical optimization. For a survey on 
complementarity problems see [4] . 

There are clearly alternative ways to solve this problem without using such 
a solver. See for recent examples [5], [T]. However nonlinear problems with 
complementary constraints occur in many circumstances as for example in a 
black oil model [3] or in problems with dissolution-precipitation. In [H] this 
latter problem was solved using a semi-smooth Newton method inspired by [S] . 
Therefore we believe that it is worthwhile to build a general purpose appro- 
priate solver for nonlinear problems with complementarity constraints, and the 
problem of a liquid-gas system with dissolution of hydrogen is a very good test 
problem. 
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5 Conclusion 

We showed how to extend Henry's law to the case where the gas phase dis- 
appears by plugging it into a phase diagram formulation. This results into the 
formulation of a two-phase (liquid-gas) two-component (water-hydrogen) system 
with possible dissolution of hydrogen as a system of nonlinear partial differential 
equations with a set of nonlinear complementarity constraints. 
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